Search for: All records

Abstract. Global Navigation Satellite Systems (GNSS) radio occultation (RO) is one of the most vital remote sensing techniques globally and of major importance for numerical weather prediction (NWP) and climate science. However, retrieving profiles of atmospheric quantities such as temperature or humidity from GNSS observations is not straightforward and dedicated algorithms still have their limitations. One of these limitations is the need for external meteorological data in the retrieval process. Various new RO missions have led to an enormous increase in data amounts and with over 10000 globally-distributed, daily profiles, RO can be considered big data nowadays. In this study, we make use of this fact by developing a new retrieval method based on a deep learning model, which only needs RO-specific quantities as an input to produce atmospheric profiles. The model is trained on almost a full year of data from COSMIC-2 and Spire RO missions, using vertical profiles of bending angle (BA) and other RO parameters as input features and operational results from a standard retrieval algorithm as target values for supervised learning. Initial results from both internal and external validation using reanalysis and radiosonde data suggest that this method produces results with an accuracy comparable to standard algorithms, while mitigating the need for external information in the retrieval process itself. These initial results serve as a starting point for further development of data-driven models for RO, which could significantly enhance the quality of RO products utilized in, e.g., climate sciences by mitigating external biases and increasing independence from other techniques. 

Abstract. The Radio Occultation Modeling EXperiment (ROMEX) is an international collaboration to test the impact of varying numbers of radio occultation (RO) profiles in operational numerical weather prediction (NWP) models. An average of 35,000 RO profiles per day for September–November 2022 from 13 different missions are being used in experiments at major NWP centers. This paper evaluates properties of ROMEX data, with emphasis on the three largest datasets: COSMIC-2 (C2), Spire, and Yunyao. The penetration rates (percent of profiles reaching different levels above the surface) of most of the ROMEX datasets are similar, with more than 80 % of all occultations reaching 2 km or lower and more than 50 % reaching 1 km or lower. The relative uncertainties of the C2, Spire, and Yunyao bending angles and refractivities are estimated using the three-cornered hat method. They are similar on the average in the region of overlap (45° S–45° N). Larger uncertainties occur in the tropics compared to higher latitudes below 20 km. Relatively small variations in longitude exist. The assimilation of ROMEX data caused small degradations in biases in several NWP models. We investigate biases in the observations by comparing them to each other and to models. C2 bending angles appear to be biased by about +0.1–0.15 % compared to Spire and other ROMEX data. These apparent biases, some of which are representativeness or sampling differences, are caused by the different orbits of C2 and other ROMEX missions around the non-spherical Earth and the associated varying radii of curvature (radius of a sphere that best fits the Earth’s surface curvature at a given location and orientation of the RO occultation plane and is used in the RO BA retrievals). 

We compare the random error statistics (uncertainties) of COSMIC (Constellation Observing System for Meteorology, Ionosphere and Climate, C1) and COSMIC-2 (C2) radio occultation (RO) bending angles and refractivities for the months of August 2006 and 2021 over the tropics and subtropics using the three-cornered hat method. The uncertainty profiles are similar for the two RO missions in the troposphere. However, a higher percentage of C2 profiles reach close to the surface in the moisture-rich tropics, an advantage of the higher signal-to-noise ratio (SNR) in C2. C2 uses signals from both GPS (Global Positioning System) and GLONASS Global Navigation System Satellites (GNSS). The GPS occultations show smaller uncertainties in the stratosphere and lower mesosphere (30–60 km) than the GLONASS occultations, a result of more accurate GPS clocks. Therefore, C2 (GPS) uncertainties are smaller than C1 uncertainties between 30–60 km while the C2 (GLONASS) uncertainties are larger than those of C1. The uncertainty profiles vary with latitude at all levels. We find that horizontal gradients in temperature and water vapor, and therefore refractivity, are the major cause of uncertainties in the tropopause region and troposphere through the violation of the assumption of spherical symmetry in the retrieval of bending angles and refractivity. 

Abstract Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇_HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇_HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇_HN| below 8 km. Significance StatementThese results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.